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Abstract. Several of the models for origin of the highest energy cosmic rays also 
predict significant neutrino fluxes. A common factor of the models is that they must 
provide sufficient power to supply the observed energy in the extragalactic component 
of the cosmic radiation. The assumption that a comparable amount of energy goes 
into high-energy neutrinos allows a model-independent estimate of the neutrino signal 
that may be expected. 



An important argument in favor of supernova explosions as the power source for 
galactic cosmic rays is the fact that kinetic energy of the ejecta supplies the right 
amount of power. The energy content of the cosmic radiation is 



where = dN/dE is the measured local flux of cosmic rays corrected for the 

effect of solar modulation. In the source region the average energy density in cosmic 
rays is related to the average production rate per unit volume, q{E), by 



where Tesc is the characteristic residence time of cosmic rays in the source region 
(for example the disk of the galaxy). The characteristic time, Tesc{E), decreases 
with energy [1] so that the observed spectrum is somewhat steeper than the source 
spectrum. 

Given an estimate of t^sc and the rate of supernova explosions, it is possible 
to estimate the fraction of energy of supernova explosions needed to maintain the 
galactic cosmic rays in steady state, assuming that supernovae provide the power, 
q{E). Assuming a rate of three supernovae per century with a kinetic energy of 10^^ 
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I INTRODUCTION 





Pe = q{E) X t^sc{E), 



(2) 



ergs per supernova, the conclusion is that an efficiency of ~ 10% for conversion of 
kinetic energy of supernova ejecta into relativistic cosmic rays would suffice. In the 

years since this coincidence was pointed out [2], a theory of cosmic-ray modified 
shocks with a high efficiency for particle acceleration has been developed [3]. A 
compelling feature of this theory is that it produces a spectral index that fits in 
well with what is observed after energy dependence of propagation is accounted 
for. 

Here we want to apply a similar analysis of energetics to cosmic rays of ex- 
tragalactic origin The reason to focus on the high energy end of the cosmic ray 
spectrum in connection with high energy neutrino astronomy is that several mod- 
els have been suggested as sources of ultra-high energy cosmic rays which would 
also be likely sources of high energy neutrinos. These include some models of active 
galactic nuclei (AGN) and some models of gamma-ray burst (GRB) sources. 

In the case of galactic cosmic rays, if the production rate of neutrinos is pro- 
portional to that of cosmic rays, qu{E) — fx q{E), we would expect a flux of 
neutrinos related to the cosmic-ray flux by 

Vd 

OC X / X 0cr-, (3) 

c 

where / is the efficiency with which the cosmic rays interact to produce neutrinos 
(either in the source region or in the interstellar medium) and Vd is the velocity 
of diffusion of cosmic rays out of the galaxy. Vd ~ Rgaiaxy I 'Tesc depends on the 
model of cosmic-ray propagation in the galaxy, but in any case Vd <C c. Thus the 
galactic neutrino flux should be suppressed by a large factor relative to the parent 
cosmic rays because of their straight-line propagation out of the galaxy. For the 
same reason, neutrinos, like gamma-rays, would be expected to show the structure 
of the galactic disk. 

For a cosmological distribution of sources of ultra-high energy cosmic rays we 
assume initially, for the rough estimates discussed below, that 



D / 3000 Mpc 
Vd ~ Rh th ~ —77, ^ c. (4) 



This approximation neglects evolution and assumes that the intergalactic magnetic 
fields are weak enough so that charged particles can reach us from distant sources in 
less than the age of the universe. The situation for super-GZK particles, where the 
maximum distance of propagation is limited to i? <IC Rh by photo-pion production, 
needs a separate discussion. The possibility of large inter-galactic magnetic fields 
and local concentrations of sources [4,5] would also have to be considered for a 
full treatment of the cosmic-ray spectrum. We use the approximation of Eq. 4 
primarily to derive an estimate of the diffuse flux of high energy neutrinos that 
may be associated with the sources of extragalactic cosmic rays. 

The outline of the paper is as follows. First we discuss the transition from galactic 
to extragalactic cosmic radiation in order to deflne an extragalactic component. 



Then we compare the power needed for the extragalactic cosmic radiation with 
that available from various potential sources. We conclude with a review of the 
predictions for detection of high energy neutrinos. 

II EXTRAGALACTIC COMPONENT OF COSMIC 

RAYS 

It is generally believed that sources of the highest energy cosmic rays are extra- 
galactic, or at least not confined to the plane of the galaxy. Indeed, there is some 
evidence for a transition from one particle population to another somewhere above 
10^^ eV. There is a trend from heavy toward fighter composition in the measure- 
ments of Xmax versus energy (see Fig. 1) and a suggestion of a hardening of the 
spectrum between 10^^ and 10^*^ eV (see Fig. 2). The picture is not as clear as it 
first appeared in the original stereo Fly's Eye result [6]. For example, the coincident 
measurements of the prototype HiRes Fly's Eye with the MIA ground array [7] , as 
shown by the gray, filled circles in Fig. 2, have a steeper slope and trend toward 
the proton curve more quickly than the original Fly's Eye data (black triangles). 
This would indicate a stronger change of composition at a somewhat lower energy. 

Another relevant piece of information is the anisotropy measurements from 
AGASA [13] and Fly's Eye [14], which can be interpreted as an enhancement of 
particles from near the direction of the galactic center in the energy bin around 
10^* eV [15]. The anisotropy disappears at higher energy. 

Although the experimental picture is still not entirely clear, in order to esti- 
mate the power we assume in what follows that the transition to cosmic rays of 
extragalactic origin occurs between 10^^ and 10^^ eV. 

Ill POWER FOR EXTRAGALACTIC COSMIC RAYS 

Measurements of the cosmic-ray spectrum above 10^^ eV are summarized in Fig. 
2. There are at least three problems that must be dealt with to estimate the power 
required to supply the extragalactic cosmic radiation. First is where the transition 
to the extragalactic component occurs in the data. As discussed above, we assume 
that the particles with ii^ > 3 x 10^^ eV are mostly of extragalactic origin and 
normalize the extragalactic component at > lO-*^^ eV. To set the scale for the 
uncertainty in this assumption, we investigate below the consequences of increasing 
or decreasing the crossover energy by half a decade. The second problem is to decide 
how to extrapolate to lower energy where the observed spectrum is likely dominated 
by cosmic rays from inside the galaxy. This is important because most of the energy 
content is fikely to be in the lower energy particles. Finally, we must deal separately 
with the super-GZK particles whose sources must be relatively nearby. We deal 
with these points in turn. 
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FIGURE 1. Average depth of shower maximum (Xmax) vs. energy [6-9] compared to calcu- 
lated [10] protons (upper curves) and iron primaries (lower curves) in two models [11,12]. 



A Normalization at high energy 

Normalizing an E'"^ differential spectrum, as is generally assumed for the sources 
of the highest energy cosmic rays [22,23], to a measured spectrum with a differential 
spectral index close to a = 3 is a highly uncertain operation. The lower heavy line 
in Fig. 2 shows a spectrum with a = 2 and an exponential cutoff at 5 x 10^^ eV 
to represent the GZK effect. The excess of data above the curve for E < 10^^ eV 
may be attributed to the high energy tail of the galactic cosmic-ray spectrum. 
Another possibility could be that we live inside a cosmologically local concentration 
of extragalactic sources of cosmic rays which propagate with an energy- dependent 
Tesc analogously to galactic cosmic rays. 

Integrating the energy content under this curve as in Eq. 1 gives for the energy 
density in cosmic rays of extragalactic origin, peg ~ 2 x 10^^^ erg/cm'^. Note that 
this result assumes the extragalactic spectrum extends down to ~ 1 GeV. Replacing 
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FIGURE 2. The high energy cosmic rays spectrum. See text for explanation of curves. Data 
are from Refs. [16-21,6] 
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Tesc ^ T// ~ 10^° yrs in Eq. 2 then leads to the estimate qeg ~ 10^^ erg/Mpc^/s. 
Shifting the normahzation point lower (or higher) by half a decade in energy would 
increase (decrease) this estimate by roughly a factor of two. This is comparable to 
the systematic differences among the different measurements of the spectrum. 

Table 1 shows what this power requirement would imply for various classes of 
potential sources. In each case, these are comparable to observed luminosities. 
Therefore all are plausible potential sources provided a mechanism exists to achieve 
Emax ~ 10^° eV. (For clusters of galaxies see Refs. [24,25]; for active galaxies [26] 
and for GRB [27,28].) 



B Extrapolation to low energy 



If the spectral index of the extragalactic source spectrum is steeper than a = 2, 
then the power requirement will be greater. Both AGN and (especially) GRB 
involve relativistic shocks. Acceleration at relativistic shocks typically produces a 
spectral index a ~ 2.25 [29-31]. The difference is particularly important in the case 
of GRB where the bulk Lorentz factor F ~ 300. Such a spectrum, with the same 
normalization at 10^^ eV is shown by the upper heavy line in Fig. 2. In general 



TABLE 1. Power per source for extragalactic cosmic 
rays 





Source density (rate) 


Power per source 


3 X 10"^ galaxies/Mpc^ 
3 X 10-*^ clusters/Mpc^ 
10-'^ AGN/Mpc^ 
1000 GRB/yr 


3 X 10^^ erg/s/galaxy 
3 X 10^^ erg/s/clustcr 
Iq44 erg/s/Active galaxy 
3 X 10^2 erg/GRB 



q{a = 2) \ E^in J (« - 2) \n{E^ax / Emin) ' 

where the normaUzation is fixed at 10^^ eV. For a ^ 2.3 and Emm ~ 1 GeV, the 
power requirement is a factor of ~ 100 greater than for a = 2.0. 

Vietri [32] argues, however, that in the case of a relativistic shock accelerating 
particles from swept up material, Emin ~ x rUp, which corresponds to Emm ^ 
100 TeV for GRB with F ~ 300. In this case the enhancement factor in Eq. 5 
is only about a factor of 3, and the power requirement increases from 3 x 10^^ to 
10^3 erg/GRB. 



C Super-GZK particles 

One possibility is that the particles above the GZK cutoff may be due to a lo- 
cal concentration of the same type of sources that produce a universal component 
[33,34]. The hne labelled Super-GZK in Fig. 2 represents a possible contribution 
from such nearby sources, assuming a = 2. The energy integral for this compo- 
nent is approximately a factor of four lower than for the corresponding universal 
contribution, but the power density to supply it is higher by a factor of five or 
ten because the distance from which sources can contribute is limited by energy 
loss due to photo-pion production. A factor of 5-10 local overdensity of sources 
of ultra-high energy cosmic rays is difficult to reconcile with other considerations. 
For example, Ref. [35] points out that most of the GRB rate comes from sources 
with z ^ 1. (See Ref. [36] for a different GRB scenario.) A local overdensity of 
more than a factor ~ 2 would also be difficult to reconcile with data on the large 
scale distribution of matter in our vicinity of the universe [37]. Given the limited 
statistics, however, this explanation cannot be ruled out if the source spectrum is 
hard enough (a < 2 [37]). 

Another possibility is that the super-GZK particles have a different origin al- 
together, being products of parton cascades generated by decay or annihilation of 
GUT-scale objects, such as topological defects [38] or massive relic particles [39]. In 
both cases the cascade consists of hadronization of partons at extremely high mass 



scale. Thus the ratio of photons and neutrinos to protons at production is large be- 
cause pions dominate the hadronic cascade. However, in the case of a cosmological 

distribution of sources, the photons will initiate electromagnetic cascades during 
propagation in the microwave background, reducing their contribution to the ob- 
served super-GZK events [40]. In the model of decaying massive relics, however, the 
predominant contribution comes from relatively nearby - particles concentrated in 
the halo of the galaxy. In this case, electromagnetic cascading will be negligible, 
and the highest energy events would, for the most part, have to be photon-initiated 
showers. It has been known for some time that the big Fly's Eye event [41], for 
which the profile is measured, looks more like a shower initiated by a proton or 
nucleus than a photon [42]. Recently, an analysis of horizontal air showers [43] 
showed that most of the observed events above the GZK cutoff cannot be initiated 
by photons. 

IV EXPECTED NEUTRINO FLUXES 

A standard technique to search for high energy neutrinos of astrophysical origin 
is to look for upward-moving muons induced by v^j^ that have penetrated the Earth. 
The signal is the convolution 

Signal ~ Area ® Na ® (Tv® (f)u, (6) 

where i?^ is the muon range in g/cm^ and Na is Avogadro's number. The range and 
cross section both increase linearly with energy into the TeV region, after which the 
rate of increase slows. Neutrinos with E^, < 100 TeV are not strongly attenuated 
by the Earth, and much of the solid angle away from the nadir remains accessible 
up to 1 PeV [44]. Thus the optimum range for z/^-induced upward muons is from 
a TeV to a PeV. Also in this energy range the muon energy loss is greater than 
minimum ionizing, which is a potential way to discriminate against the background 
of atmospheric neutrinos, which have a steeply falling spectrum. In what follows, 
I use the cross sections from Ref. [44], taking account of absorption by the Earth, 
to estimate z/^-induced signals in a kilometer-scale detector. A minimum path- 
length of 0.5 km inside the detector is required, which corresponds to a threshold 
of approximately > 100 GeV in water. 

A Generic estimate 

Using the normalization of § III-A and assuming a ~ 2.0 for the neutrinos as well 
as the cosmic rays, one estimates a signal of / x 30 events/km^/yr [45], where / is 
the efficiency for production of neutrinos relative to cosmic rays. For f — 1 this is 
essentially the Waxman-Bahcall upper bound [22], which applies for sources that 
are transparent to neutrons. If evolution is included, this estimate may be increased 
by factor of about five if sources are assumed to evolve similarly to the rate of star 



formation [22,46]. This is because the ultra-high energy protons from high redshift 
would be attenuated by photoproduction and pair production while the neutrinos 

would not. Thus the neutrino flux would be greater for a given cosmic-ray flux at 
the normalization point. In addition, if the spectrum of extragalactic cosmic rays 
to which the normalization is made has a > 2.0, then the corresponding estimated 
neutrino flux in the TeV to PeV range and the corresponding signal could also be 
larger. 



Rates of neutrino-induced upward muons in a km detector 
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FIGURE 3. Summary of predicted neutrino signals in a kilometer-scale detector. The diffuse 
signal per year per steradian per logarithmic energy interval is plotted as a function of E^,. 
Numbers indicate events per year above atmospheric background (see text). 



B Predictions of models 



Fig. 3 shows the rates of neutrino- induced upward muons predicted in various 
specific models. These are diffuse fluxes summed over 27r sr with absorption in 
the Earth accounted for. The background induced by atmospheric is shown 
as the heavy sohd hne. The flattening of the atmospheric background at 1 PeV 
is due to an assumed prompt component of neutrinos from charm decay [47] , the 
level of which is rather uncertain [48]. The numbers indicate events per year above 
atmospheric background, which I have estimated by taking the integral of the signal 
flux above the energy where it crosses the atmospheric background. This crossover 
is generally in the range of 100 - 300 TeV. Present upper limits from Prejus [49], 
Baikal [50] and AMANDA [51] rule out one of the original AGN core models [52], 
and the model of Ref. [53] (shown as AGN core in Fig. 3) is marginally allowed [49] 
or marginally ruled out [51]. The GRB curves in Fig. 3 are calculated from Eq. 6 
starting from the unevolved GRB neutrino flux plotted in Ref. [22]. I assumed 
source evolution from Ref. [22] to obtain the GRB-evolved neutrino spectrum. 

The GRB model [54], has an interesting feature, which illustrates the important 
point that the neutrino spectrum need not have the same shape as the parent 
spectrum of accelerated protons [55]. In this case the target for pion production 
is assumed to be the X-ray/7-ray photons in the expanding flreball, which have 
approximately an E~'^ differential spectrum above a characteristic energy e^, and 
an E^^ spectrum below the break. For protons with energy sufficiently high to 
photoproduce on photons with < ei,, the resulting neutrino spectrum follows the 
parent proton spectrum. For lower energy protons, the density of target photons 
above threshold for photo-pion production decreases as the proton energy decreases. 
The result of this convolution is a break (steepening) in the neutrino spectrum, 
which, for the parameters of Ref. [54] occurs at Ei, ~ 100 TeV. Since both and 
a^, are increasing with energy, this is where the signal peaks. 

The key to detecting a diffuse signal above the atmospheric background is to look 
for an excess of events of high energy. Typically, muons with energies of several 
TeV and higher will radiate one or two bursts per kilometer of water in which they 
deposit some 10% of their energy. Such bursting tracks may be a useful signature. 

C Point sources 

With point sources, the atmospheric background can be reduced to an extent 
that depends on the angular resolution of the detector. As an example. Fig. 4 
shows the signal that would be generated by a neutrino flux normalized to the level 
of the TeV 7-ray emission during the extended high-state of Mrk501 in 1997 [56-58], 
assuming (py^+u^ = (p-y, (speciflcaUy using the flt of Ref. [56]). The atmospheric 
background is calculated assuming that the detector has a pointing resolution of 
1° space angle, and the angle between the detected muon and the neutrino that 
produced it has been accounted for. The integral of the signal shown in the flgure 



would give ~ 30 events per year above atmospheric background. Unfortunately, 
MrkSOl is not normally such an intense source [59] as it was for half of 1997. 
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FIGURE 4. Predicted signal of z/-induced muons from a source with = ^-y{Mvk 501, high). 
(See text.) 

It is interesting that there is at least one AGN-blazar model in which a neutrino 
flux comparable to the 7-ray flux would be expected [60]. In this model the high 
energy photons and neutrinos come from decay of pions produced when blobs of 
ultra-relativistic gas collide with the interstellar medium of the host galaxy near 
the central engine of the AGN. Farther out the material slows down and merges 
into the jets that extend to large distance. 



D Ultra-high energy neutrinos 

Some of the "top-down" models for the ultra-high energy cosmic rays [61-63] 
predict an intensity of neutrinos with energies in the PeV range and higher that 
is somewhat above the neutrinos photo-produced by ultra-high energy cosmic rays 
on the microwave background (GZK neutrinos). The rates of i/^u-induced upward 
neutrinos shown for those models in Fig. 3 are strongly attenuated by absorption 
in the Earth (although regenerated could emerge with degraded energy [65]). 
However, the energy is high enough so that charged- current interactions of z/g from 
above could be identified by the large electromagnetic cascades they would produce 



inside a kilometer-scale detector [66]. For example, the rate of downward events 
predicted by the model from Ref. [63] would be ~ 10 per year per km^ even though 
the integrated rate of upward i/^-induced muons shown in Fig. 3 is ^ 1. 

E Concluding comment 

Estimates of the type discussed here - either based on energetics and a generic 
relation to extra-galactic cosmic rays, or on specific models - set the scale for a 
high energy neutrino telescope. The estimates are of the order of tens of events per 
year in a kilometer-scale detector. While the hope is that new kinds of "hidden" 
neutrino sources [67] will be discovered, the detector should be designed to find 
signals from likely known sources at this level. 
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